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ABSTRACT 

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula 
of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a 
hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin 
of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code 
cmfgen, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be 
highly peculiar (Ry ~4.7), and obtain a surprisingly high luminosity log(L/L Q ) = 6.5 ± 0.2, corresponding to a present-day mass of ~ 150M G . 
The high effective temperature of 52.2 ± 2.5kK might be explained by chemically homogeneous evolution - suggested to be the key process in 
the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are 
unprecedented for classical Wolf-Rayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the 
origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or 
(ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date. 
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1. Introduction 

Recent claims for the existence of very massive stars (VMS) 
of up to 300M o at the centre of young clusters like R136 
( |Crowther et al.||2010| l seem to link the formation of such ob- 
jects to environments in the centres of massive clusters with ~ 
10 4 - 10 s M Q . This finding thus appears to support massive star- 
formation models that invoke competitive accretion and possibly 
even merging in dense clusters (e.g. |Bonnell et al.|2004| . Such 
scenarios emerged after it became clear that the disk-accretion 
scenario - commonly applied to low- and intermediate star for- 
mation - had problems explaining the formation of stars with 
masses above 10M Q , as radiation pressure on dust grains would 
halt and reverse the accretion flow onto the central object (e.g. 
Yorke & Krue gel|1977[ ). However, recent multi -dimensional hy- 
drodynamical monolithic collapse simulations indicate that mas- 
sive stars may form via disk accretion after all (e.g. Kuiper et al. 
|2010| l. This illustrates that the discussion on clustered vs. iso- 
lated massive star formation is still completely open (see de Wit 
[eTaTI < |2005> ; |Parker & Goodwin| ( |2007T >; |Lamb et al.| p010) ; 
Bressert et al. in prep.). 

In this Letter, we present evidence that VFTS 682 
(RA05h38m55.51sDEC-69°04'26.72") in the Large 
Magellanic Cloud (LMC), found in isolation from any 
nearby massive, visible cluster, is one of the most massive stars 
known. The object suffers from high dust obscuration and is 
located in, or in the line-of-sight towards, an active star-forming 
region. On the basis of a mid-infrared (mid-IR) excess, Gruendl 
& Chu (2009) classified the star as a probable young stellar 




* Based on observations at the European Southern Observatory Very 
Large Telescope in programme 182.D-0222. 



Fig. 1. Combined YJK S image from the VISTA Magellanic 
Clouds survey ( [Cioni et al.|201 l| l. The arrows correspond to 5 pc 
in the Northern and Eastern directions. The projected distance of 
VFTS 682 to the cluster R136 is 29 pc (for an assumed LMC dis- 
tance of 50kpc). 



object (YSO). Th e VLT-FLAMES Tarantula Survey (VFTS) 
( |Evans et al.|201 1| ) identified it as a new Wolf-Rayet (WR) star, 
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Table 1. Stellar parameters of VFTS 682 



Parameter 


Value 


7"* 


54.5 ± 3 kK 


Teff 


52.2 ± 2.5 kK 


log(L/L ) 


6.5 ±0.2 


logCM/VZ/Moyr- 1 ) 


-4.13 ±0.2 


P 


1.6 


f 


0.25 




2600 ± 200 km/s 


Y 1 


0.45 ±0.1 


M v 


-6.83 ±0.12 mag 



Notes. (1) Y is the helium mass fraction 

at a projected distance approximately 30 pc northeast of R136. 
Its spectrum was classified as WN5h and is similar to those of 
the very luminous stars in the core of R 136, with strong H, He n, 
and N iv emission lines. In this Letter, we present a photometric 
and spectroscopic analysis of VFTS 682 to investigate its origin. 
The star is relatively faint in the optical, which we argue is the 
result of significant reddening (with Av ~4.5), implying a high 
intrinsic luminosity and a mass of order 150 M Q . The sheer 
presence of such a massive star outside R136 (and apparently 
isolated from any notable cluster) poses the question of whether 
it was ejected from R 136 or if it was formed in isolation instead. 

2. Spectroscopic Analysis of VFTS 682 

The analysis is based on spectroscopic observations (/14000 - 
7000) from MEDUSA mode of VLT-FLAMES ( |Evans et al.] 
|201 1). We fir st compared observations taken at different epochs 
fEvans et al.||20ll) as to identify potential shifts in the radial 
velocity (RV) as a result of binarity. To this purpose, we used 
the Niv/14060 emission line. With detection probabilities of 
96% for a period P<10d, 76% for 10d<P<lyr, and 28% for 
lyr<P<5yr, we can basically exclude that VFTS 682 is a short 
period binary. We found the Nv/14944 line to be particularly 
useful for RV determinations, as this line remained stable for 
different wind velocity laws. With particular emphasis on the 
centroid of the line we obtained 300 + 10 km/s. If we give more 
weight to the red wing, we obtain 315 + 15 km/s, with a similar 
number for the He n /J4686 line. 

For the spectral analysis we use the non-LTE model atmo- 
sphere code cmfgen (Hillier & Miller 1998 ) which has been de- 
veloped to provide accurate abundances, stellar parameters, and 
ionising fluxes for stars dominated by wind emission lines. We 
use atomic models containing H i, He 1-n, C m-iv, N m-v, O m-vi, 
Si rv, P iv- v, S iv- vi, Fe iv-vii and Ni iv-vi. We made the following 
assumptions. We adopt a metallicity of half solar with respect to 
Asplund et al. ( 2005 i, and assume N and He to have been en- 
riched by the CNO process. We adopt a 12-fold enhancement 
of the N mass fraction. We employ a B parameter of 1 .6 for the 
wind velocity law and a wind volume filling factor of / = 0.25, 
assuming that the clumping starts at 10 km/s. In order to estimate 
the error bars in the stellar temperature (T+) and mass-loss rate 
(M), we computed a grid of models around the preferred values. 

The core temperature (T*) at high optical depth of 54500 ± 
3000 K is based on the Niv ,14060 and N v doublet /14604 and 
4620 in emission (Fig.|2]i. Niv is also sensitive to the mass-loss 
rate, which is obtained from other lines. The hydrogen and the 
helium lines are used to determine both the helium abundance 
(45% by mass) and the mass-loss rate, with log(M/M Q yr~ l ) = 
-4.4 + 0.2, where the volume filling factor corrected mass- 
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Fig. 2. Relative flux vs. wavelength in A. Grey solid line: 
MEDUSA spectrum of VFTS 682. Black dashed line: model 
spectrum. 
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Fig. 3. Log flux vs. log wavelength. Squares: B, V, J, H, K s 
photometry, Diamonds: spitzer SAGE. Circles: spitzer 5.8 ymi 
and 8.0 /jm from |Gruendl & Chu|p009] l. Solid line: SED with 
anomalous Ry parameter. Dashed line: Standard Ry parameter. 
Dashed-Dotted line: attempt to reproduce the NIR excess with a 
1500 K black body. 

loss rate is given in Table [T] This value is larger than the 
\og(M I M Q yr~ l ) = -4.69 obtained from the mass-loss recipe of 
|Vink et aL (2001 1 for non-clumped winds. The values can be 
brought into agreement for a volume filling factor of / = 0. 1 . In 
the absence of UV spectroscopy of the /U548, /il551 Civ dou- 
blet, we estimate the terminal velocity ( — 2600 km/s) from 
the broadening of the Hen and H a lines. An overview of the 
stellar parameters and abundances is provided in Table [T] 

The stellar parameters of VFTS 682 are similar to those of 
R136a3 (with a log(L/L G ) = 6.58) ( |Crowther et aLl[20T0> . We 
note that there are other luminous WNh objects in 30 Dor, for 
instance, VFTS 617 (=Br 88) which has the same spectral type 
as VFTS 682 but with weaker emission lines. The spectral anal- 
ysis of the additional VFTS WNh stars is currently in progress 
(Bestenlehner et al. in prep.). 



3. Luminosity and Reddening 

In its spectral appearance, VFTS 682 is almost identical to that 
of R136a3 in the core of the dense cluster R136. The luminous 
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Color-Magnitude Diagram 
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Fig. 4. VFTS objects. Rectangles are WN-stars. Triangles are 
WC-stars. The arrows in the top right corner are for two different 
Ry values, with a length corresponding to Ay = 1 . The longer 
lines on the left indicate the transformations between reddened 
and unreddened CMD positions. 



objects in this cluster are thought to be the most massive stars 
known, with initial masses up to 300 M Q ( Crowther et al.|2010| l. 
As VFTS 682 is apparently isolated, and shows no signs of bi- 
narity, it offers the opportunity to study such a luminous WN5h 
object in isolation. The precise luminosity of VFTS 682 is thus 
of key relevance for this work. Its derivation is however compli- 
cated by the large extinction. 

Matching the optical photometry with a standard extinction 
law, i.e., a 'Galactic average' extinction parameter Ry = 3.1, 
we obtain a luminosity of log(L/L ) = 5.7 + 0.2. With this rela- 
tively low luminosity we can explain the optical flux distribution, 
but we fail to match the observed spectral energy distribution 
(SED) at longer wavelengths. In the near-mid infrared (IR) the 
observed flux is ~ 3 times higher than the corresponding model 
flux. Below, we show that we can successfully explain this 'ex- 
cess' IR emission, including its exact shape, with a reddening 
parameter Ry = 4.7, which leads to a much higher stellar lumi- 
nosity. This effect is visualised in the colour-magnitude diagram 
(CMD) of Fig. [4] We note that high Ry parameters, that devi- 
ate from the Galactic average, are not extra-ordinary for massive 
stars in the LMC (e.g. |Bonanos et al.|2011| l 



3.1. Modelling the near IR flux 

Because the extinction in the IR is very low, we can avoid the 
problems resulting from the uncertain Ry parameter, by choos- 
ing the near-IR photometry as our flux standard. In more detail, 
we derive the luminosity of VFTS 682 by matching B, and V 
from Parker| ([1993 ), and K s from the "InfraRed Survey Facility 
(IRSF ) Magellanic Clouds Point Source Catalog" ( |Kato et ai] 
2007| > (see Tab.[2]l. For this purpose we extract the intrinsic B, V, 



and K s magnitudes from our model, using appropriate filter func- 
tions. Based on the resulting values of E(B - V) = 0.94 ± 0.02, 
and E(V - K s ) = 3.9, we derive R v = A V /E(B - V) for two 
oft-used extinction laws. 

To determine R v , we use the relations Ry = 1.1994 x E(V - 
K S )/E(B — V) — 0.183, as inferred from the extinction law by 
ICardelli et al.] ( [T9"89l ), and Ry = 1.12xE(V-K s )/E(B-V)+0.016, 



from Fitzpatrick] ( [1999] ). This way, we obtain Ry = 4.7 ± 0.1, 
Ay = 4.45 ± 0.12, and A Ks = 0.55 + 0.15, where the uncer- 
tainties refer to the differences in the adopted extinction laws. 
As A v is much larger than A# s , the derived luminosity mostly 
relies on K s , while the observed values of B, and V mostly de- 
termine Ry. Small uncertainties in B - V will thus mainly affect 
Rv, in a way that the absolute values of B, and V stay the same. 
As Ry = Ay/E(B - V), this introduces an anti-correlation be- 
tween the derived Ry, and E(B - V), where Ay is preserved. The 
resulting luminosity is log(L/L Q ) = 6.5 ± 0.2, where the error 
includes uncertainties in the stellar temperature, mass-loss rate, 
and extinction. The resulting SED fit is presented in Fig. [3] 

In an independent analysis of the available UBVIJHK pho- 
tometry, using the CHORIZOS code ( |Mafz-Apellaniz| |2007| l 
for CHi-square cOde for parameteRized modeling and char- 
acterization of photometry and Spectrophotometry, we obtain 
E(B - V) = 0.98 ± 0.03, and R v = 4.55 ± 0.17, which, due to 
the aforementioned anti-correlation between E(B - V) and Ry, 
results in an almost unchanged value of Ay = 4.46 ± 0.06. The 
resulting luminosity, and high R v are thus confirmed. 

Table [2] gives a summary of the available photometry for 
VFTS 682. In analogy to our previous findings, the differences 
between the optical photometry of |Evans et al.| ( |201 l| l, and 
Parker ( 1993 1 mainly affect Ry, but not the derived luminos- 
ity. Based on the WFI photometry by Evans et al. ( 201 1) , we 
obtain Ry = 5.2 ± 0.1 with E(B - V) = 0.84 ± 0.02. The 
near-IR photometry however has a direct influence on our re- 
sults. Based on the fainter Two Micron All Sky Survey (2MASS) 
photometry the derived luminosity decreases by about 0.1 dex to 
log(L/L ) = 6.4. In this work, we use the IRSF values, as they 
connect better to the observed spitzer MIR photometry. 

If we alternatively adopt the standard reddening parameter 
Ry = 3.1, and attribute the NIR excess to a second component, 
such as a cool star and/or a warm (1000-2000 K) dust component 
(see Fig. [3j, we obtain a luminosity of only log(L) = 5.7 + 0.2. 
We note that just one extra component (cool star and/or dust 
black body) cannot fit the entire SED slope, which implies that 
we would need a multitude of additional components. These 
would need to conspire in such a way as to precisely follow 
our SED. Although we cannot disprove such a configuration, it 
would be highly contrived. 

To summarise, the most likely luminosity of VFTS 682 is 
log(L) = 6.5 + 0.2, which would support the high luminosities 
that have been derived by Crowther et al. ( 2010| l for the VMS in 
the core of R 136. 



3.2. The MIR excess 

We also investigate the issue of whether VFTS 682 has a mid- 
IR excess in spitzer data. We compared two sources of MIR 
photometry, the "Surveying the Agents of a Galaxy's Evolution 
(SAGE) IRAC Catalog" for point sources (|Meixner et al.|2006] 



06) 
)9». 

09) 



as well as the recent YSO catalog of |Gruendl & Chu|( |2009) 
In contrast to the SAGE IRAC CatalogT lGruendl & Chu| ( |2009) 
included spatially more extended sources. The MIR photometry 
data from both sources are listed in Table [2] Gruendl & Chul 
(2009) detected a mid-IR excess at 5.8 and 8.0 ymi and cate- 
gorised VFTS 682 as a YSO candidate. We are not able to con- 
firm the MIR excess on the basis of the SAGE IRAC Catalog 
photometry, but the difference in the two data sources may be 
explained in the following way. 

The MIR point-source data of Meixner et al. (2006) basically 
follow the slope of our optical-to-NIR fit continued into the MIR, 
whilst the larger scale MIR data of Gru endl & Chu| ( f2"009| ) show a 
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Table 2. Apparent Magnitude of VFTS 682 
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Fig. 5. OGLE III V and I lightcurves for VFTS 682 during 2001- 
2009. 



clear MIR excess. Whereas the MIR excess could be consistent 
with an extended circumstellar shell, such as recently detected 
in spitzer 24 micron images of the Galactic Centre WNh star 
WR102ka (Barniske et al. 2008), there is no evidence for such 
a shell (or bowshock) in spitzer 8 micron images of VFTS 682 
at similar angular scales as for WR102ka. However, given the 
resolution limit of ~ 2 arcsec the IR excess might be due to an 
unresolved circumstellar shell with a diameter of < 0.5 pc. 

4. Variability on timescales of years 

To study the stability of VFTS 682 we show Optical 
Gravitational Lensing Experiment (OGLE-III) lightcurves 
( |Udalski et aT][2008l ) in V and I in Fig. [5] On top of the short- 
term jitter, the object clearly dimmed by ~0.1 mag in both the 
V and I band during the years 2006-2009. Furthermore, the ob- 
ject shows an increase in the K-band by ~0.15 mag during 2010 
(Evans et al. in prep.). This kind of long-term variability is un- 
precedented for Wolf-Rayet stars, and is more characteristic for 
Luminous Blue Variables (LBVs). However, we note that the 
nature of these changes is not the same as those of bona-fide 
LBVs with S Doradus cycles of 1-2 mag variations (Humphreys 
|& Davidson|i994] i. 



5. Discussion on the origin of VFTS 682 

VFTS 682 is located in an active star-formation region of 30 Dor 
( Johansson et al.fl998 i. It is not anywhere close to the core of 
R136, or - as far as we can tell - to any other nearby cluster. This 
might be considered a surprise as VMS are normally found in the 



cores of large clusters, such as the Arches cluster or R136. By 
contrast, VFTS 682 is rather isolated at a projected distance of 

29 pc from R136. This raises the question of whether the object 
formed in situ or is a slow runaway object from R136 instead. 

In order to address these issues, we provide some velocity es- 
timates. The measured RV shift of VFTS 682 is v r « 300km/s, 
which is somewaht higher than the average velocity in 30 Dor 
(~ 270±10km/s, Sana et al. in prep, |Bosch et aLl[2009| . If 
VFTS 682 is a runaway from R136, it would require a tan- 
gential velocity of 30km/s to appear at a projected distance of 

30 pc within approximately 1 Myr. Together with an RV offset 
of ~30km/s, VFTS 682 would then require a true velocity of 
~40 km/s, i.e. at the lower range of velocities for classical OB 



runaway stars (e.g. Philp et al. 



19961. Still, this would make it 



the most massive runaway star known to date. 

In case the object is a runaway star, a bow shock might po- 
tentially be visible as VFTS 682 is surrounded by dust clouds. 
Whilst there is currently no evidence for a bow shock, one of 
the possible explanations for the MIR excess is the presence of 
a bow-shock (e.g. Gvaramadze e t al.| 2010i. Alternative expla- 
nations for the MIR excess may involve its (line-of-sight) asso- 
ciation with an active star forming region, or a nebula formed 
by vigorous mass loss since the object started to burn hydro- 
gen in its core. In this context, it may or may not be relevant 
that VFTS 682 shows slow photometric variations suggestive of 
LBVs that may suffer from mass ejections. 

The most probable luminosity of VFTS 682 is log(L) = 6.5 + 
0.2. The question is what stellar mass the object corresponds 
to, and what is its most likely evolutionary age. The high tem- 
perature places the object right on the zero-age main sequence 
(ZAMS), which can be best understood as a result of chemically 
homogeneous evolution (CHE). Such an evolution has also been 
proposed for two WR stars in the SMC (Marti ns et""aTp 009). 
Note that this type of evolution has been suggested for the pro- 
duction of long gamma-ray bursts ( |Yoon & La nger 20051. In 
order to obtain a meaningful mass limit, we employ the recent 
mass-luminosity relationships by Grafener et al. ( 201 1 j > for ho- 
mogeneous hydrogen burners. Utilising the derived He abun- 
dance (Y=0.45), the most likely present-day mass is ~150M G . 
In detailed stellar evolution models of |Brott et ah] ( |201 \\ and 
Friedrich et al. (in prep), CHE is achieved through rapid rotation, 
with v™' > 200km/s. A value of Y=0.45 is obtained at an age 
of 1-1.4 Myrs and the initial mass would be of ~ 120 - 21OM , 
where the mass range is defined through the uncertainty in the 
luminosity. 

Finally, from our analysis it is clear that VFTS 682 is a very 
massive object. It is often assumed that such massive objects 
can only be formed in dense cluster environments, where they 
are normally found. The apparent isolation of VFTS 682 may 
thus represent an interesting challenge for dynamical ejection 
scenarios and/or massive star formation theory. 
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